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Neutrino- induced production (neutrinoproduction) of photons and pions from nucleons and nuclei 
is important for the interpretation of neutrino-osciUation experiments, as these photons and pions 
are potential backgrounds in the MiniBooNE experiment [A. A. Aquilar-Arevalo et al. (MiniBooNE 
Collaboration), Phys. Rev. Lett. 100, 032301 (2008)]. These processes are studied at intermediate 
energies, where the A (1232) resonance becomes important. The Lorentz-covariant effective field 
theory, which is the framework used in this series of studies, contains nucleons, pions. As, isoscalar 
scalar (a) and vector (lj) fields, and isovector vector (p) fields. The Lagrangian exhibits a nonlinear 
reahzation of (approximate) SU{2)l ® SU(2)b. chiral symmetry and incorporates vector meson 
dominance. In this paper, we focus on setting up the framework. Power counting for vertices 
and Feynman diagrams is explained. Because of the buih-in symmetries, the vector current is 
automatically conserved, and the axial-vector current is partially conserved. To calibrate the axial- 
vector transition current {N -ir^ A), pion production from the nucleon is used as a benchmark and 
compared to bubble-chamber data from Argonne and Brookhaven National Laboratories. At low 
energies, the convergence of our power-counting scheme is investigated, and next-to-leading-order 
tree-level corrections are found to be small. 
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PACS numbers: 25.30.Pt; 24.10.Jv; ll.30.Rd; 12.15.Ji 



I. INTRODUCTION 

Neutrinoproduction of photons and pions from nucle- 
ons and nuclei plays an important role in the interpre- 
tation of neutrino-oscillation experiments, such as Mini- 
BooNE The neutral current (NC) 7r° and photon 
production produce detector signals that resemble those 
of the desired signals. Currently, it is still a question 
whether NC photon production might explain the excess 
events seen at low reconstructed neutrino energies in the 
MiniBooNE experiment, which the MicroBooNE experi- 
ment plans to answer [2|. Moreover, pion absorption after 
production could lead to events that mimic quasielastic 
scattering. 

Ultimately, the calculations must be done on nuclei, 
which are the primary detector materials in oscillation 
experiments. To separate the many-body effects from 
the reaction mechanism and to calibrate the elementary 
amplitude, we study charged current (CC) and NC pion 
production from free nucleons in this work, which serves 
as the benchmark. Moreover, NC photon production, 
which is not a topic under intense investigation, is stud- 
ied within this calibrated framework. In future papers, 
we will include the electroweak response of the nuclear 
many-body system to discuss the productions from nuclei 
in the same framework. 

Here we use a recently proposed Lorentz-covariant 



meson-baryon effective field theory (EFT) that was orig- 
inally motivated by the nuclear many-body problem '3|- 
[To| . (This formalism is often called quantum hadrody- 
namics or QHD.) This QHD EFT includes all the rel- 
evant symmetries of the underlying QCD; in particu- 
lar, the approximate, spontaneously broken SU{2)l <E) 
SU(2)ii chiral symmetry is realized nonlinearly. The mo- 
tivation for this EFT and some calculated results are dis- 
cussed in Refs. In this EFT, we have the A 
resonance consistently incorporated as an explicit degree 
of freedom, while respecting the underlying symmetries 
of QCD noted earlier. (The generation of mesons and 
the A resonance through pion-pion interactions and pion- 
nucleon interactions has been investigated in pll |22|.) 
We are concerned with the intermediate-energy region 
^^Lab ^ 0.5 GeV), where the resonant behavior of the 
A becomes important. The details about introducing A 
degree of freedom, the full Lagrangian, and electroweak 
interactions in this model have been presented in (23), [2j] . 
The well-known pathologies associated with introducing 
A are not relevant in the context of EFT. The cou- 
plings to electroweak fields are included using the ex- 
ternal field technique [1^ , which allows us to deduce the 
electroweak currents. Because of the approximate sym- 
metries in the Lagrangian, the vector currents are au- 
tomatically conserved and the axial-vector currents are 
partially conserved. Form factors are generated within 
the theory by vector meson dominance (VMD), which 
allows us to avoid introducing phenomenological form 
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factors and makes current conservation manifest. -"^ We 
discuss the power counting of both vertices and diagrams 
on and off resonance and consistently keep ah tree-level 
diagrams through next-to-leading order. Explicit power- 
counting of loop diagrams in this EFT has been discussed 
in Refs. |17Hl9l |. Here the contributions of the loops are 
assumed to be (mostly) saturated by heavy mesons and 
the A resonance, so the couplings of contact interactions 
are viewed as being renormalized. The mesons' role in 
effective field theory has also been investigated in (27l[28| . 

One major goal of this work is to calibrate electroweak 
interactions on the nucleon level. It is typically assumed 
that the vector part of the N ^ A transition current is 
well constrained by electromagnetic interactions [l^, HO] ■ 
The uncertainty is in the axial-vector part of the cur- 
rent, which is determined by fitting to Argonne National 
Laboratory (ANL) [3l| and Brookhaven National Labo- 
ratory (BNL) [s^l bubble-chamber data. The data have 
large error bars, which leads to sig nificant model depen- 
dence in the fitted results [26|, [s^, Isii Here we choose 
one recently fitted parametrization [33!| and use it to de- 
termine the constants of our VMD parametrization (but 
note that our basis of currents is different from the con- 
ventional one as used in [2^ [sl, HJ]). In addition, we 
make use of other form factors, the ones in [2^ for exam- 
ple. We then compare results of using different current 
basis and form factors with the data at low and interme- 
diate neutrino energies. 

There have been numerous earlier studies of neutrino- 
production of pions from nucleons in the resonance region 
[H, [H, HO, m ii-lii- Theybasically fall into two cate- 
gories. In the first one [2^, HO, IH, HI, [H, |4l| resonance 
dominance above intermediate energy is assumed. The 
contributions of resonances are summed incoherently and 
hence it is difficult to determine the interference effect. In 
the second category, [1^ HI, H^, El] , the contributions 
are summed coherently including the background, since 
either an effective Hamiltonian or Lagrangian is utilized. 

Our approach belongs to the second category, while dif- 
ferences from other models should be mentioned. First, 
there exists a finite energy range in which EFT is valid, 
so we insist on low-energy calculations. However, a differ- 
ent attitude has been taken, for example, in Refs. (26ll43 |. 
in which the Born approximation based on an effective 
Lagrangian has been extrapolated to the region of sev- 
eral GeV. Second, we have discussed the consequence of 
higher-order contact terms. ^ Naturally, these contribu- 
tions should obey naive power counting [4J,'45'|; however, 
some of them may play an important role in scattering 
from nuclei. Third, electroweak interactions of nucle- 
ons are calibrated in this work while the strong interac- 
tion has been calibrated to nuclear properties. This is 



a unique feature that is absent in other models target- 
ing the production from free nucleons only. Furthermore, 
the calibration on the nucleon's electroweak interaction 
impacts the strong interaction. For example, the p tt tt 
coupling, introduced because of VMD in the pion's vec- 
tor current, gives rise to an interesting contribution in 
the two-body axial current in a many-body calculation 
[46| . In our theory with A, it can be quite interesting 
to investigate similar consequences, for example, the A's 
role in the two-body current, in which meson-dominance 
couplings can give rise to relevant interactions. 

This article is organized as follows: in Sec. Ulland IIIIl 
we introduce our Lagrangian without and with A, and 
we calculate several current matrix elements that will be 
useful for the subsequent Feynman diagram calculations. 
The theory involving A is emphasized. Then the tran- 
sition current basis and form factors are discussed care- 
fully. In Sec. IIV| we discuss our calculations for the CC 
and NC pion production and for the NC photon produc- 
tion. After that, we show our results in Sec.|Vl Whenever 
possible, we compare our results with available data and 
present our analysis. Finally, our conclusions are sum- 
marized in Sec. I VII 

In the Appendixes, we present the necessary informa- 
tion about chiral symmetry and electroweak interactions 
in QHD EFT, form factor calculations, power counting 
for the diagram with A, and kinematics. 



II. LAGRANGIAN WITHOUT A (1232) 

In this work, the metric g^i, = diag(l, — 1, — 1, — l)^i,. 
The convention for the Levi-Civita symbol e^""^ is 
g0i23 _ -[^ -y^g have introduced upper and lower isospin 
indices [1^ H^]. In this section, we focus on the La- 
grangian without A and study various matrix elements: 
{N\V^,Al,J^\N) and (A^; 7r|y;, 4, J^f |iV). Definitions 
of fields and currents can be found in Appendix [K\ 



A. Power counting and the Lagrangian 

The organization of interaction terms is based on 
power counting EE] and naive dimensional anal- 

ysis (NDA) [4J, |45| . We associate with each interaction 
term an index: i) = d + n/2 + b. Here d is the number of 
derivatives (small momentum transfer) in the interaction, 
n is the number of fermion fields, and b is the number of 
heavy-meson fields. The Lagrangian is well developed in 
Refs. [13, mill 113 ■ We begin with the Lagrangian 
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^ Meson dominance generates form factors for contact pion- 
production vertices automatically, as shown in diagram (f) in 
Fig. [l] In other approaches, for example [26j , these form factors 
are introduced by hand, which requires specific relations between 
the nucleon vector current and the pion vector current form fac- 
tors. This is explained in Sees . [ITB] and ITVl 

Some of these terms have also been discussed in Ref. 42]; how- 
ever, the interpretation of these terms is different here from that 




^N{0^3) = N{ij^[d^ + igpPf, + tg,V^] + gAl^^l"^ a,, ~ M + gs<p)N 



r 



» ■ J c J • t:^ iTJ-ci Z a — i • ^ d\uPi^] + *5/> : Pi' ■ The superscripts and denote 

dn IS defined m Eq. jA6|, = - [d^ - tv^ + iy(s)„), ^i'^ • • at I • i ■ ■ 

^ , ^^ n , 1 i \r o TA o TA J the isospm. Next is a purely mesonic piece: 

and the held tensors are Vfj,i, = Of_iV^ ~ Oi, and = 



>C„(,^4) = -d^q^d'^ct>+^f''rr[d,,U{d'^U)^] + -f^mlTr{U + U^ -2) 

-\MP,.Pn - \ V^'^V,, + ^ (^Tr(F(+)^V^.) + i fl!'^V,}j . (2) 



We only show the kinematic terms and photon couplings 
to the vector fields. The latter are used to generate VMD. 
Other ly — 3 and v = 4 terms in £meson(e^4) are impor- 
tant for describing the bulk properties of nuclear many- 
body systems and can be found in [2^ [U, |4^, |49| . 



The only manifest chiral-symmetry breaking is through 
the nonzero pion mass. Other chiral-symmetry violating 
terms and multiple pion interactions are not considered 
in this calculation. Finally, we have 



jONM->=i) = ^ Njp{2f3^°^d,fr + P^^^d^pWt^'' + p^^^-i^d^F^-^nN - c^i Tr(^^^+) + TT{a^d,F^-^n 
+ ^3 Tr {a^^i 

+ ^Nra'^Nj,^^+'^NrNTy{a'^p^,) + '-f^Nra-^NV,, . (3) 

I 

Note that Cn,-jt{0=4) is not a complete list of all possible they are the only relevant z> = 4 terms for NC photon pro- 

;) = 4 interaction terms. The terms listed in the first duction. Further discussion will be given in Sees. IIV CI 

two rows generate the form factors of currents for nucle- and IV Dl 
ons and pious, gp-^^ is used for VMD. Special attention 
should be given to the ci, ei, cip, and eip couplings, since 
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B. Contributions to current matrix elements from 
irreducible diagrams 

To calculate various current matrix elements, we need 
to understand the background fields in terms of elec- 
troweak boson fields; this connection is given in Ap- 
pendix |Al Based on the Lagrangian, we can calculate the 
matrix elements (A^|V;j,Aj^, Jf^\N) and (iV;7r|V;i, A^, 
J^\N) [diagram (f) in Fig. [1] at tree level; loops are 
not included; only diagrams with contact structure are 

^ The expressions for the currents Usted below differ from those in 
Refs. [l(ll46l | because contributions from non-minimal and vector 
meson-dominance terms are included here. 

I 



We can also use this procedure to expand the axial- 
vector current in powers of using the Lagrangian con- 
stants gA and Z?^^'*- In fact, we can improve on this by 
including the axial-vector meson (ai^) contribution to 
the matrix elements, which would arise from the inter- 



included'^. Because of VMD, we can extrapolate the cur- 
rent to nonzero [13, [IS] ■ The results are given below, 
and the explicit calculations are shown in Appendix |B] 
Note that q'' is defined as the incoming momentum trans- 
fer at the vertex; in terms of initial and final nucleon mo- 
menta, q^' = p'^^f - p'^i- Similarly, q^' -fp^^ = pj'^^ + k!^ for 
pion production. 

First, the matrix elements of the nucleon's vector and 
baryon current, and the axial- vector current in pion pro- 
duction are the following: 



(4) 
(5) 



(6) 



(7) 
(8) 



(9) 
(10) 

I 

actions: ga^N^f^^^aif^N and Ca, Tr {F}~^^"'alf,„) . Here 
oi/i = aitf,T'/2 and ai^^ = d^ai^, - d^ai^, where a^^ 
are the fields of the ai meson (with its mass denoted as 
TOai — 1.26 GeV). Then we obtain 



{N,B\J^\N,A)=5iuf[^^ + 25Fl 



{N,B-t:,j,KW^\N,A) - (^b\^\A)ujYu, 



u., = SgUfTBi^iq)ui , 



2M 

V.ryid,,.. 1. ^2\i- {i- k-K)g^LU - [q- K)^,qi, 



g^, + 25F^-'^''{{q~K 



'-f-{B\-\A)u} 



(Jfj.,yiq 

2M 



2A(i) + 25F^^'^\{q - k^f) 



g ■ (g - fcx) 
[q - KY 



-y- {B\ — \A)ufrA-^t^{q,k^)ui . 

Jit ^ 



Here nip = 0.776 GeV, 0.783 GeV, bF = F{q^) ~ F{0) (also true for other form factors), and 



f; 



V,md 



F, 



V.md 



S,md 



— q — 



2A( 



1) fp9p q^ 



9-y q^ - ml 



\ , =-1.35, ^ = 2.48 
/ 9t 



g^ g2 _ ^2 



= ; 1 



/?(") , 25. 



3.97 q^ — ITT, 



, A(i) = 1.85, /„ = 3.04 , 



j) , = -1.40, ^ = 3.95 , 

iJ 9i 



pS,md ^ U 2;^(0) _ ^!9v9^\ ^ ^(0) ^ _Q ^ _o^9 



3^7 q^ - ml 
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{N,B\A^^\N,A) 
{N,B,^,j\V;\N,A) 



7^Uj = {B\ — \A)ufTAt,{q)ui , 



(11) 



{B\ — \A) UfTvT,^,{q,K)u^ , 
,(1) 2ca,ga,q'^ 



GT{q') ^ 9A- P'a' in - > 9A - 1.26, /J^^^ = 2.27, c,,g,, = 3.85 



(12) 
(13) 



For the pion's vector current form factor Q, 



{2K - g)^ + 25Fr{q') [ k,^ - ^ 



Qp-niT q 

57 9^ - "l2 



57 



= 1.20 . 



i^\Pvfi{q,kn 



(14) 



To determine the couphngs in Eqs. 0, ©, ©, (ITUl) . 
([13]) and (fl^. we compare our results with the fitted 
form factors [131 • We require that the behavior of our 
vector- and baryon-meson-dominance form factors near 
= be close to that of the fitted form factors [50| . 
The nucleon's axial-vector current used to fit our G'^'^ is 
parametrized as GA{q^) = .9^/(1 — q^ /Ma)"^ with gA = 
1.26 and Ma = 1.05 GeV [51]. As shown in Ref. I^, 
the form factors due to vector meson dominance become 
inadequate at « 0.3 GeV^. This is also true for the 
axial-vector's parametrization. This indicates that the 
EFT Lagrangian is only applicable for Ei ^ 0.5 GeV in 
lepton-nucleon interactions, above which exceeds the 
limit. This will be clarified in the kinematical analysis of 
Sec. TO 



III. LAGRANGIAN INVOLVING A(1232) 
A. Lagrangian 

Two remarks are in order here [2^[23]: First, the the- 
ory is self-consistent with general interactions involving 
tp^] second, the so-called off-shell couplings, which have 
the form ^f^ip'^, d^j^ip^ ^tp^ and d^j^ip^ , can be consid- 
ered as redundant. For the chiral symmetry realization, 
A*" belong to an / = 3/2 multiplet [a = (±3/2, ±1/2)]. 
Moreover in the power counting of vertices, the A is 
counted in the same way as nucleons. 

Consider first £a [v ^i): 



- 4^ A, p^^a^-'A' - 4^ A.V^.a^-'A' - ^ A,v,,a^-^A' + ^ A, A^ Tr(S' a,) . (15) 
4m 4m m m 

I 



This is essentially a copy of the corresponding Lagrangian To produce the -s-> A transition currents, we con- 
fer nucleons. struct the following Lagrangians (j> ^ 4): 
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C^,N,. = hAA'"'T^'^ai^NA + C.C. , 

•^A,JV,bg 



(16) 



M 
+ C.C. 



M 



M2 



+ C.C. . 



(17) 
(18) 



Here Tj 



jl which are (complex conju- 



gate of) Clebsch-Gordan coefficients. 



Based on the Lagrangians, we find (noting that 



B. Transition currents 

We can express the transition current's matrix element 
as follows: 

{A,a,pA\V"'{A^n\N,A,pN) = 

T^'^UAaipA)r^'lA)il)^NiPN) . (19) 



r7 



2ciA(g') 



-flA 9 



g2 - ml 



M2 



M 



(20) 



(,«< - g'^^q') - ^ (g«7'' - 5"'' ^ - ^ ^"a'^'^i?. , (21) 



M2 



CiA(9^ 



ClA 



Qa + 



CiAp 



2^^ q'^-mj ' 



i = 1,3,6, 



1.21 , C3A = -0.61 , C6A = -0.078 , 



ClAp 
57 



-4.58 



C3Ap 
57 



2.32 



C6Ap 
9-y 



0.30 . 



(22) 



Similar to the CiA{q ), we can introduce axial- vector me- higher-order terms contained in 5hA{q ) = hA{q ) — hA 

son exchange into the axial transition current, which conserve the ax;ial transition current. With the lim- 

leads to a structure for the diA{q^) similar to that of the ited information about manifest chiral-symmetry break- 

CiA{q^)- There is one subtlety associated with the rcaliza- ing, we ignore this subtlety and still use the form of the 

tion of hA{q^)'- with our Lagrangian, we have a pion-pole Cia(9^) to parametrize hA{q^)- 
contribution associated with the Ha coupling, and all the 
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/ia(q^) =hA+ llAai 



g2 _ 



, hA = 1.40 , /lAai = -3.98 



(23) 



, i = 2,4,7, 



9 '""ai 

d2A = -0.087, d4A - 0.20, d7A = -0.04 , d2Aai = 0.25 , 



rf4Aai = -0.58 , d7Aai = 0.12 



(24) 



To determine the coefficients in the transition form fac- (A, ^|jcc+l^' [— ^^2/3 Ua(pA) (r^^ + F^'') u(pAr)] 
tors shown in Eqs. (22) and we compare ours is 

with one of the conventional form factors used in the ht- 
erature. In Refs. [l^, [s^ for example, the definition for 



Ma (pa) 



M 



C\ 



cl 



(5"'^ ^ - q"r ) + ^il-PA 5"^ - 9>^) + C^g'^' + % q^'q'^ 



M 



u{pn) 



(25) 



We use the "Adler parametrization" 35] in Ref. [3^ to fit 
our meson-dominance form factors. Now supposing the 
baryons are on shell, we can represent the conventional 
basis as linear combinations of our basis. For example, 

q^a'^'iq^j^ = (m - Af)(g"7^ - 5"'' ^j-f'^ 

- {q>^A - 1 ■ PAg'-ni' - {qy^, - q ■ PNg'^^i^ ■ (26) 

A similar relation holds with 7^ deleted on both sides and 
{rn~M) changed to (m + A/). We can obtain the relation 
between form factors associated with the two bases: 



ClA 



C3A 



m-M {CX + CX) 



2M 



(27) 



C6A 



16 



(28) 



(29) 





(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



FIG. 1: Feynman diagrams for pion production. Here, C 
stands for various types of currents including vector, axial- 
vector, and baryon currents. Some diagrams may be zero 
for some specific type of current. For example, diagrams (a) 
and (b) will not contribute for the (isoscalar) baryon current. 
Diagram (e) will be zero for the axial- vector current. The 
pion-pole contributions to the axial current in diagrams (a) 
(b) (c) (d) and (f) are included in the vertex functions of the 
currents. 



diA 



m - 



MCf 



2M 



, d7A = \ 7: 



3 
2 



(30) 



We assume that these relations hold away from the res- 
onance. It can be shown that, at low energy, the dif- 
ferences in observables due to using the two bases, with 
these relations applied, are negligible. Moreover, the q^ 
dependence of these c^a a-nd diA form factors can be real- 
ized in terms of meson dominance. We then require that 



the meson-dominance form factors be as close as possible 
to the ones indicated in Eqs. (|27l) to ((30l) . and we get the 
couplings shown in Eqs. (1^^ . (1^51) . and ([M)) . However, 
we should expect the leading-order meson-dominance ex- 
pressions would fail above w 0.3 GeV^. 
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IV. FEYNMAN DIAGRAMS 

Tree-level Feynman diagrams for pion production due 
to the vector current, axial- vector current, and baryon 
current are shown in Fig. [T] In this section, we calculate 
different matrix elements for pion production and photon 
production. The Feynman diagrams for photon produc- 
tion can be viewed as diagrams in Fig.[l]with an outgoing 
TT line changed to a 7 line. It turns out that diagram (e) 
in Fig. [1] is negligible in NC photon production, since it 
is associated with 1 — 4 sin^ 9w • 

First let us outline the calculation of the interaction 
amplitude M . Consider CC pion production (in the one- 
weak-boson-exchange approximation) : 

M ^ 4V2Gj.Kd(^S^)(^r"''')- ■ (31) 

where i = +1, —1. In Eq. (pT|) . Gp is the Fermi con- 
stant, Vud is the CKM matrix element corresponding to u 
and d quark mixing, and {J^^°''^'^^^) = {NB, i:j\J'^^\NA) 
(The definitions of currents can be found in Appendix [XI) 

{J^l^^) = {liJ)\J Liti\vi{i'i)) is the well-known leptonic- 
charged-current matrix element. For NC pion pro- 
duction, we need to set Vud = 1, ( j|''°'^^*'')7r — > 
{J^^fn^. and {J^^) ^ in Eq. Here 

('^n'cI} leptonic-neutral-current matrix element, 

and {J^^c'^'^^)^ = {NB,'Kj\J%^\NA). For NC photon 



production, we have an expression similar to that of 
NC pion production with ( J^'^''^^)^^ ^> {Jp^cf''^)-!, while 

Now consider the power counting for ( j('"^'')'')jr(-y) in 
Eq. pip for various processes. The order of the diagram 
{v) is counted as ^ = 2L + 2- En/2 + J2^#iih-2), 
where L is the number of loops. En is the number of ex- 
ternal baryon lines, Vi = di + ni/2 + bi is the order of the 
vertex (j)) mentioned in Sec. Ill A[ and #i is the number 
of times that particular vertex appears. However, there is 
a subtlety related with power counting of diagrams with 
A, which has been carefully discussed in Ref. [52]. Com- 
pared to the normal power counting mentioned above, in 
which the baryon propagator scales as 1/0{Q), for dia- 
grams involving one A in the s channel, we take i> — !■ 1 
in the resonance regime and i> + 1 away from the 
resonance. Details can be found in Appendix [Cl 

Finally, conservation of vector current, conservation of 
baryon current, and partial conservation of axial- vector 
current can be easily checked for the matrix elements 
shown in the following. 



A. Diagram (a) and (b) 

Diagram (a) and (b) in Fig. [T] lead to currents: 



I 



where fc^r is the outgoing pion's momentum. 
Ty^^ji^^{jp]q,Pi) are defined in Eq. (jl9p while A's momen- 
tum IS p = q + Pi. ry"^-j{pf;q,p) = 7°r|;'(^^(p; -q,p/)7° 
while p = —q + Pf. In the following, we always have this 
definition of F. The A's propagator, SFfj.u{p), is shown 
in Appendix |D] The subscript j denotes the isospin 
of the outgoing pion. For vector current, in diagram 
(a) Vnr !^ 3 in the lower-energy region and i/j, 5^ 1 in 
the resonance region; in diagram (b) Vnr 1^ 3. For 
axial-vector current, in diagram (a) Vnr 1^ 2, ^ 0; 
in diagram (b) i>nr ^2. In the power counting, the 



higher-order terms in v come from including form factors 
at the vertices. Moreover, the baryon current matrix 
element is zero {{J^)Tr = 0) in both diagrams. 

Now we examine the NC matrix element {Jj^cf'^'^)'/- 
First, based on the relations given in Eq. (|A18|) . we define 

^%^{p;q,Pi)^ Q -sin^6'^^ r^^(p;g,K) ir^^(p;(7,pO : 

(33) 

Then we find 



i^Nc)! = eTlgT^'^^ Uf el{k)Ty°'{pf;~k,p)SFo,f3ip)T'lj'ip;q,Pr)ut 

+eT°'gTl^UfT^j^ {pf]q,p)SFaf}{p)^v^{p]-k,p^) el{k)ui , (34) 
I 

where k is the outgoing photon's momentum and e^(fc) diagram (a) Vnr ^ 4, i^j. ^ 2; in diagram (b) Vnr ^ 4. For 
is its polarization. For the vector current in the NC, in 
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the axial- vector current, in diagram (a) Vnr 3, t'r 5^ 1; 
in diagram (b) z/„r ^ 3- 



B. Diagrams (c) and (d) 

These two diagrams lead to currents: 



(V^^iA^n). = --^ {B\''-^UA)uf M'SApK^^^iq^H - ^ {B\-^\A}ufr^y^^^{q)SF{p) M'u.. • (35) 



^gA 



2 2 



2 2 



For the nucleon propagator, p = q+Pi in diagram (c) and For the baryon current we just need to change ^Ty{q) 

p = -q + Pf in diagram (d). r^^^^(g) has been defined to ^^{q) in Eq. (|35|) . and i/ ^ 1. 

in Eq. (|4]). For both currents in both diagrams ^ 1. For NC photon production, we get 



Ncn 



c 



) + ^r^(-fc)j5Fb) 



A r 



C 



I -Sin' 9^) r^(g) + ir^(<z) 



-euf 



^-sinH^)T^y{q) + ^r'Xiq) 



-^-Isin'e^r'i.iq) 



xSF{p)el{k) 



where we use the shorthand {^)b — (^1^1^)- 
three currents, power counting gives v ^ 1. However, this 
naive power counting does not give an accurate compari- 
son between the A contributions and the N contributions 
at low energies, as we discuss later. 



J 



C. Diagrams (e) and (f) 



The two diagrams lead to a vector current 



(36) 



■U}{i- f^^h'u, + ^ {B\^\A)ufr^^iq,k^)u, 



(37) 



Here, Py{q, fc^) is defined in Eq. (|14p . Ty^{q, fc^) is de- For the axial- vector current, diagram (e) does not con- 

fined in Eq. (|12p . and v ~^ \. tribute, and we find 
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life 



{B\-\A)ufV^^^{q, K)u, + {B\-\A) 



f 



g2 - ml 2M 



SI 



1 



-IKl 1 



4 M2 



q2 _ 



(38) 



Here, r^^(g,fc^) is given in Eq. dH). The terms in the A^' B" + A" B^' . 

first row lead to ^ 1 contributions. The contributions For the baryon current, diagrams (c) and (f) do not 

due to K^, and ki are at = 2. We use values fitted contribute: {Jg)-K ~ 0. 

in [5^ for these couplings. In the last row, A^f^B'^^ — For the NC photon production matrix clement we find 



Here v = 3; for j/ < 3, there are no contact vertices 
contributing in this channel. By power counting, we ex- 
pect that at low energy, these terms can be neglected 
compared to the v = 1 terms. However according to 
Ref. |42l] . these terms may play an important role in co- 
herent photon production. Meanwhile, it is claimed in 
Ref. [42] that the origin of these contact vertices is the 
anomalous interactions of the uj and p. But they can also 
be induced by the off-shell terms in the A Lagrangian. 
Moreover, we can construct meson-dominance terms by 
using the interaction terms in the last row of Eq. ([3]) and 
photon- meson coupling in Eq. which leads to differ- 
ent off-shell behavior of the vertex compared to that of 
the anomaly term. 




FIG. 2: The configuration in the isobaric frame. 



A. Kinematics 



V. RESULTS 



In this section, after introducing the kinematics, we 
discuss our results for CC and NC pion production, and 
also NC photon production, and compare them with 
available data whenever possible. 



Fig. [2] shows the configuration in the isobaric frame, 
i.e., the center-of-mass frame of the final nucleon and 
pion. The momenta are measured in this frame, except 
those labeled as p^, which are measured in the Lab frame 
with the initial nucleon being static. Detailed analysis of 
the kinematics is given in Appendix [El The expression 
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for the total cross section is 
1 



IMI- 



32Af„ (27r)5 K + E,,f Ip^f 



IMI- 



1 



^ ^ ■; — r-, — r dflT^dM^„dQ^ , 

64M2 (2^)5 E^ + E,,f ^ ' 

(39) 



where |Mp is the average of total interaction ampli- 
tude squared. Based on the equations in Appendix [El 
we can make the following estimates. For CC pion 
production, when E^ = 0.4 (0.5) GeV, (M^„)„aa; ~ 
1.17 (1.24) GeV,g^<,^ « 0.2 (0.3) GeV^ We can see 
that above E^ = 0.4 GeV, the interaction begins to 
be dominated by the A resonance. However, when 
= 0.75 GeV, {M^n)max ~ 1.4 GeV, and higher res- 
onances, for example Pii(1440), may play a role. The 
exception is + p — > ii~ + p + -it~^: only 7 = 3/2 
can contribute, and the next resonance in this channel 
is the A(1600), which is accessible only when i?^ ^ 
1.8 GeV. For NC pion production and photon production 
{Ei; ^ 0.2 GeV), when E^; = 0.3 (0.5) GeV, (M^nUa. ~ 
1.2 (1.35) GeV, Q^„^ ~ o.l (0.3) GeV^ Here above = 
0.3 GeV, the interaction begins to be dominated by the 
A. However, when Ej; = 0.6 GeV, {M^n)max ~ 1.4 GeV, 
and higher resonances may play a role. 

From this analysis, we expect our EFT to be valid at 
E^f ^ 0.5 GeV, since only the A resonance can be ex- 
cited, and Q2 ^ 0.3 GeV^ where meson dominance works 
for various currents' form factors [2^. To go beyond 
this energy regime when we show our results, we require 
Mttu ^ 1.4 GeV and use phenomenological form factors 
that work when ^ 0.3 GeV^ 



plied. In these calculations, F™'^, G"^'^, ca, and are 
substituted by the form factors in the literature. The re- 
sults of our framework, i.e. using the meson-dominance 
form factor born out of the Lagrangian, are shown as 
"MDFF" calculations, and these are extrapolated beyond 
0.5 GeV limit also. The extrapolations of both CFF and 
MDFF calculations enable us first to compare our result 
with similar calculations in (26| , ^ and second to see how 
meson-dominance form factors fail at higher energy. By 
comparing CFF with MDFF calculations, we can see in 
the MDFF calculation that the meson-dominance form 
factors are inadequate for reproducing the conventional 
form factors above E^, = 0.5 GeV (although it seems that 
MDFF results are closer to the data). Hence in the 
following Fig. [51 we only show the MDFF results with 
E^ s$ 0.5 ^ 0.6 GeV, for which M^„ ^ 1.4 GeV holds 
automatically. Since we believe the EFT is applicable in 
this low-energy regime, in these plots, we show results in- 
cluding Feynman diagrams up to order v = 1 and v = 2. 



CFF only A ' 
HFFonly A 
MDFF only A 
CFF upto v=1 
HFF upto v=1 
-MDFF upto v=1 

ANL with cut I 
BNL without cut 



o 
o 




0.8 1 
Ev(GeV) 



B. CC pion production 



In this section, we compare calculated cross sections 



of CC pion production with ANL 31] and BNL (32| 
measurements. In both experiments, the targets are hy- 
drogen and deuterium. [All the other experiments use 
much heavier nuclear targets in (anti) neutrino scatter- 
ing, and to explain this, we must examine many-body 
effects.] The beam is composed of muon neutrinos, the 
average energy of which is 1 and 1.6 GeV for ANL and 
BNL respectively. In the ANL data, there is a cut on 
the invariant mass of the pion and final nucleon sys- 
tem: Mttk ^ 1.4 GeV; no such cut is applied in the BNL 
data. Based on the previous phase-space analysis, this 
cut clearly reduces the number of events when E^, is above 
'-^0.5 GeV. This can be seen by comparing the two data 
sets in three different channels shown in Figs. [3l and [H 
the ANL data lie systemically below the BNL data. Since 
the data stretch above 0.5 GeV, in Figs. [3l and [4l we show 
the "CFF" results [using the conventional form factor in 
p]] and the "HFF" results [using the form factor in ^ 
with the reduced C^(0)], with the M,™ constraint ap- 



FIG. 3: (Color online). Total cross section for ly^ +p — 5- 
H~ +P + TT^. "Only A" indicates that only diagrams with 
A (both s and u channels) are included. The "up to u = 1" 
category includes all the diagrams at leading order. The CFF 
calculations are done with one of the conventional form factors 
33| . The HFF calculations make use of form factor used in 
3 with the reduced C^(0). The MDFF calculations are 
based on the EFT Lagrangian with meson dominance. In the 
ANL data, Mnn ^ 1.4 GeV is applied, while no such cut is 
applied in the BNL data. For all calculations, Af^„ ^ 1.4 
GeV is applied. 



In Fig. [21 we show the data and calculations for 
+ p — > /i^ + p + 7r+. As mentioned above, in the 
"CFF only A" calculation, we make use of one set of 
conventional form factors and include the Feynman dia- 
grams with the A in both s and u channels. In the "CFF 



* The calculation in [Sal without reduction of C^{0) should be 
close to the CFF calculation [33j |. although the details of the 
form factors are different. 
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up to u = 1" calculation, we use the same form factors 
and include all the Feynman diagrams up to leading or- 
der. These two calculations are quite similar to those 
done in Ref. without reducing C^. Indeed, our re- 
sults are consistent with theirs. (In Ref. ^261, only the s 
channel contribution is included in the calculation with 
"only A.") Next, we show two different HFF calcula- 
tions: one with only A (in the s and u channels) and 
the other with all the diagrams up to = 1. Finally, we 
also show two MDFF calculations up to different order, 
so that we can compare the MDFF approach with the 
CFF approach. 

First, we can see that both CFF and MDFF calcula- 
tions with only A diagrams are consistent with the data 
at ^ 0.5 GeV. Introducing other diagrams up to order 

= 1 is still allowed by the data at low energy, although 
they indeed increase the cross section noticeably. Second, 
in Ref. [26|, a reduced C^(0) is introduced, primarily to 
reduce the calculated cross sections above — 1 GeV, 
which can be seen by comparing CFF calculations with 
HFF calculations. However, since we are only concerned 
with the El, ^ 0.5 GeV region, in which we see satisfac- 
tory agreement between our calculations and the data, 
we will keep the C^(0) fitted from the A's free width. 
Furthermore, in the original spectrum-averaged da /dQ^ 
data of ANL [3l|, the contributions from E^ < 0.5 GeV 
neutrinos are excluded, so comparing calculations with 
data at low energy is not feasible at this stage, and we 
will not show our da/dQ^ here. 

In Fig. m we show the data and calculations for 

+ n — > /i^ + n + 7r+ and + n — > /i^ + p + tt*^. 
We can see that the situations in these two processes are 
quite similar to that in Fig. [3l the results of the CFF 
and MDFF approaches are consistent with the data at 
low energy. Again the differences between the two ap- 
proaches with the same diagrams begin to show up when 
the neutrino energy goes beyond 0.5 GeV. Although the 
pion production is still dominated by the A, if we com- 
pare cross sections (from the same calculation) in Figs. [3] 
and m we see that other diagrams introduce significant 
contributions and violate the naive estimate of the ra- 
tio of the three channels' cross sections based on isospin 
symmetry and A dominance. Moreover, the reduction 
of significantly reduces the cross section in these two 
channels if we compare the two HFF calculations with 
the corresponding CFF calculations. 

In Fig. [5l we begin to investigate the convergence of 
our calculations in different channels in neutrino and an- 
tineutrino scattering. We show the MDFF calculations 
based on our EFT Lagrangian up to different orders. We 
see that the power counting makes sense systematically 
in different channels: including N intermediate state and 
contact terms up to = 1 changes the "only A" calcula- 
tion non-negligibly. Far below resonance, the A contri- 
bution is less important compared to that in other dia- 
grams, and it begins to dominate around 0.4 GeV. This is 
consistent with the power counting discussed in Sec. IIVI 
Moreover, the v = 2 terms do not change the "up to 



3 F CFF only A' 
HFF only A 
MDFFonlyA 
2-5 - CFFuptov=1 
HFF up to v=1 
J- 2 JVIDFF up to v=1 
E ANL with cut 




0.2 0.4 0.6 0.8 1 1.2 1.4 



Ev(GeV) 



CFF only A' ' 

HFFonlyA 

- MDFFonlyA 

CFFuptov=1 

- HFFuptov=1 

MDFFuptov=1 

ANL with cut 1 — 1 — 1 


1 


1 


1 








BNL without cut 1 — X — i ' 
- vn->u~Dii° .. 1— , 










1 j<<iM*^ 1 ' 1 1 


1 1 



0.2 0.4 0.6 0.8 1 1.2 1.4 



Ev(GeV) 

FIG. 4: (Color online) Total cross section for (a) -I- n — Y 
-\- n + TT^ and (b) v^, + n — > fi~ -I- p -I- tt". In the ANL 
data, Mttu ^ 1.4 GeV is applied, while no such cut is applied 
in the BNL data. The curves are defined as in Fig. [3] 

1/ — 1" results significantly. All the calculations of neu- 
trino scattering are consistent with the limited data from 
ANL. We can see that the cross section for antineutrino 
scattering is generally smaller than that of neutrino scat- 
tering, due to the relative sign chosen between V^'^ and 
in the Feynman diagrams having A. The sign be- 
tween V'^ and A''^ in other diagrams is well defined in 
our Lagrangian. The relative sign between A's contribu- 
tion and other diagrams' is also well determined by the 
relation between Ha and in Eq. ([2^ . although it has 
been investigated phenomenologically in Ref. [26j. 



C. NC pion production 

In this section, we discuss the results for NC pion pro- 
duction in (anti)neutrino scattering. In Fig. [SJ the re- 
sults in the MDFF approach are shown for calculations 
including diagrams of different orders. The channels are 
explained in each plot. Since all of the available data 
for NC pion production are spectrum-averaged, and neu- 
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FIG. 5: (Color online) Total cross section for CC pion production due to neutrino and antineutrino scattering off nucleons. 
Here "only A" indicates that only diagrams with A (both s and u channels) are included, "up Xo v = 1" includes all the 
diagrams at leading order, and "up to v = 2" includes higher-order contact terms, whose couplings are from Ref. [s^. In the 
ANL data, M-^n ^ 1.4 GeV. For calculations, M-^n ^ 1-4 GeV is applied. 



trinos with i?^ ^ 0.5 GeV have a small weight in such 
analyses, we do not compare our results with data. Here 
we focus on the convergence of our calculations; intro- 
ducing the V = 2 terms does not change the total cross 
section significantly. However, we also see the viola- 
tion of isospin symmetry in the "up to v = 1" and 
"up to V — 2" calculations in each plot, if we com- 
pare each pair of channels in Fig. [6l In principle, if 
there is no baryon current contribution in NC produc- 



tion, we should see that the two channels in each plot 
yield the same results. For example, isospin symme- 
try implies (p,7rO|yO'^,^OA'|p) = (n, 7r°|y°'^, ^O'^jn) and 
(P,7r0|j^|p) = -(n,7rO|J^|n). So with "only A," we 
can not see the difference between the two cross sec- 
tions, since the (isoscalar) baryon current cannot induce 
transitions from iV to A. After introducing nonresonant 
diagrams, wc would expect them to be different, as con- 
firmed in the first plot in Fig.[5]for example. This analysis 
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FIG. 6: (Color online). Total cross section for NC tt production due to neutrino and antineutrino scattering off nucleons. The 
curves are defined as in Fig. [S] and the channels are also indicated. 



can be applied to other channels, and we clearly see the 
nonresonant contributions. 



D. NC photon production 

In this section "we focus on NC photon production. The 
results are sho'wn in Fig. [71 Besides NC 7r° production, 
this process is another important background in neutrino 
experiments. One important difference bet"ween NC pho- 
ton production and CC and NC pion production, is that 
all of the v = 2 terms do not contribute in this pro- 
cess. Therefore, we include the t-wo = 3 terms in NC 
photon production, namely the ei and ci couplings in 
Eq. besides terms due to the form factors. More- 

over, these t-wo couplings are singled out in Ref. [43| as 
the lo"w-energy manifestations of anomalous interactions 
involving p and uj, and they are believed to give impor- 
tant contributions in coherent photon production from 
nuclei. Here "we also investigate the consequences of these 
t'wo couplings. We emphasize that from the EFT per- 
spective, the only 'way to determine these t"wo couplings 
is by comparing the final theoretical result "with data, 
rather than by calculating them from anomalous inter- 
actions, "which are not necessarily the only high-energy 



physics contributing to these two operators. For exam- 
ple, as we discussed before, an off-shell coupling between 
N, TT, and A can introduce the same matrix element as 
the ci term. Changing the off-shell couplings would also 
change the contact term to make the theory independent 
of the choice of off-shell couplings. Nevertheless, to per- 
form concrete calculations without precise information on 
the coupling strengths, we use the values from Ref. 
(ci = 1.5, ei = 0.8). 

We can see the convergence of our calculations in 
Figs. [71 The two couplings introduced in the "up to 
1/ = 3" calculations increase the total cross section in 
both channels for both neutrino and antineutrino scatter- 
ing, although the change is quite small. This constructive 
behavior is consistent with the results in Ref. j-i^ ). 

Naive power counting, however, does not give an ac- 
curate comparison between the A contributions and the 
N contributions at very low energy. First, the neutron 
does not have an electric charge, so its current should 
appear at higher order than the naive estimate would in- 
dicate. Second, for the proton, due to the cancellation 
between the baryon current and the vector current, the 
neutral current is mainly composed of the axial-vector 
current, which reduces the strength of the neutral cur- 
rent. Because of these two effects, the contributions of 
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FIG. 7: (Color online) Total cross section for NC photon production due to neutrino and antineutrino scattering off nucleons. 
Here "only A" indicates that only diagrams with A (both s and u channels) are included, "up to i/ = 1" includes all the 
diagrams at leading order, and "up to ;^ = 3" includes higher-order diagrams. 



Compton-like diagrams are smaller than the power count- 
ing indicates. 

VI. SUMMARY 

Neutrinoproduction of photons and pions from nu- 
cleons and nuclei produce important backgrounds in 
neutrino-oscillation experiments and therefore must be 
understood quantitatively. In this work, we stud- 
ied the productions from free nucleons in a Lorentz- 
covariant, chirally invariant, meson-baryon EFT. For 
(anti)neutrino energy around 0.5 GeV, the A resonance 
is important. We therefore included the A degrees of 
freedom explicitly in our EFT Lagrangian, in a manner 
that is consistent with both Lorentz covariance and chiral 
symmetry. 

It is well known that in a Lagrangian with a finite num- 
ber of interaction terms, including the A as a Rarita- 
Schwinger field leads to inconsistencies for strong cou- 
plings, strong fields, or large field variations. In a mod- 
ern EFT with an infinite number of interaction terms, 
however, these pathologies can be removed, if we work 
at low energies with weak boson fields. This is clarified 
in our previous work [2^. Ambiguous and so-called off- 
shell couplings involving A have also been shown to be 
redundant in the modern EFT framework, because these 
couplings produce terms that can be absorbed into the 
contact terms in the EFT Lagrangian. Thus the A res- 
onance can be introduced in our EFT Lagrangian in a 
consistent way. 

Because of the symmetries built into our Lagrangian, 
the vector and baryon currents arc conserved and the 
axial-vector currents are partially conserved automati- 
cally, which is not true in some other approaches to this 
problem (with special constraints among different form 
factors having to be introduced by hand to conserve vec- 
tor current in other approaches). Needless to say, the 
conserved vector and baryon currents are crucial for com- 



puting photon production. We discussed in detail how 
the meson-dominance mechanism works in our matrix el- 
ement calculations, which is the key ingredient in current 
conservation. By using vector and axial- vector transition 
currents that were calibrated in pion production at high 
energies, we found results for pion production at lower 
energies (i?J/^^ ^ 0.5 GeV) that are consistent with the 
(limited) data. This is also true when vertices described 
by meson dominance were used. On the other hand, the 
couplings introduced to generate meson dominance are 
relevant in other problems. For example, the interac- 
tions in Eq. (|18p lead to a proper description of the vector 
transition current at nonzero and meanwhile it is rele- 
vant to two-body currents: suppose a photon is absorbed 
by one nucleon producing a A which then interacts with 
other nucleon through the interactions mentioned above. 

Moreover, we studied the convergence of our power- 
counting scheme at low energies (where the A needs 
to be counted differently in different energy regions) 
and found that next-to-leading-order tree-level correc- 
tions are small. This power counting scheme is different 
from the canonical one, because it can be used in nuclear 
many body problems. For example, the lowest order in 
this scheme is the mean-field approximation, if the cal- 
culation is done for the property of the nuclear ground 
states. The discussion on this can be found in |17Hl9j . 
It is certainly interesting to see how the power count- 
ing that we have for scattering off nucleons works in the 
scattering off nuclei. 

Finally, we computed NC photon production and ex- 
plored the power counting in this problem. The difference 
between the NC photon production and pion production 
is that, at ^ = 2, no diagrams contribute in the pho- 
ton case, while there are several in pion production. So 
we proceeded to include ly = 3 diagrams induced by two 
contact interactions, Ci and ei terms. They have been 
studied in [l^l, and they are believed to be the low en- 
ergy manifestation of anomalous p and uj interactions. 
We pointed out the existence of other sources including 
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ofF-shell couplings of A and possible meson-dominance 
terms. Nevertheless, by using two coupling strengths cal- 
ibrated to anomalous p and w interactions , we found 
that, at least for a nucleon target, their contributions are 
very small, as expected based on power counting. 

We are currently using this QHD EFT framework to 
study the electroweak response of the nuclear many-body 
system, so that we can extend our results to photon and 
pion neutrinoproduction from nuclei, which are the true 
targets in existing neutrino-oscillation experiments. 
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Appendix A: Chiral symmetry and electroweak 
interactions in QHD EFT 

Details on chiral symmetry and electroweak interac- 
tions in QHD EFT can be found in [H, [H. We intro- 
duce background fields including = y'^^Tijl (isovec- 
tor vector), m^^^s^ (isoscalar vector), and a'^ = ^^^^Tijl 
(isovector axial- vector) , where i — x,y,z or + 1,0,-1. 
They couple to the corresponding currents in QCD. 
We define r'^ = + a^, and = \/f - a^. Un- 
der SU{2)l <8) SU{2)]i (g) U{1)b symmetry transforma- 
tions, these fields should change in the following way: 
Ll^'L'< +iL^^'L'< {L = cxp[-i0u{x) ^]}, r'' ^ 
i?r^i?t + iRdt^R"! {R = exp[-iem{x) ^]}, and vj^^ ^ 



d^^O. Here, 9Li{x), 9jii{x), and 9{x) are the ro- 
tation angles. We can construct field strength tensors: 
htiu = d^l^ - dj^ ^ LfLfii^L\ and Jr^^ and 

fsiw are constructed in the same way. 

Now we discuss nonlinear transformations of dynami- 
cal degrees freedom in our model: 



U = exp 



T^i (x) ;, 



^ = Vu = exp 



iC/i?t , (Al) 
^ L^h^ = h^R'' , (A2) 



hvnh'^ — ihdnh) , 



(A3) 



-> ha^h^ , (A4) 



f(±) 



d^U- ilf^U + iUr^ L df^UR^ , (A5) 

exp He{x)B]h^ {8^,4') . (A6) 
-i\a^ , 0.1,] hv^^h'' , (A7) 
^^fL^u^±^fR^..e ^hF^^t^h^ , ( A8) 
dxFl^t^ + i[vx , Fj^t^] ^ hdxFl^^h^ (.A9) 



In the preceding equations, are the generators of re- 
ducible representations of SU{2). The « 93 MeV 
is the pion-decay constant. We generically label non- 
Goldstone isospin multiplets including the nucleon, p me- 
son, and A by ■00 = (A^A, Pi, ^a)a- B is the baryon num- 
ber of the particle. The transformations of the isospin 
and chiral singlets and are trivial. The dual field 

tensors are defined as p^^^^^ = e^'^°'^F^^\ which have 
the same chiral transformations as the ordinary field ten- 
sors. The objects shown here are the building blocks for 
constructing the Lagrangian. 

Electroweak interactions of quarks in the Standard 
Model [H, m m, [11] determine the form of the back- 
ground fields in terms of the vector bosons Wj^, Z^^ and 
A,,: 



lu. = 



.9 . 2 
■ sm I 



cos ( 

-1 T+i 



-iT-l 



9 T 



COSf 



1 , 

-e - -, — 

2 COS 0,, 



cos ( 
9 



9 . 2 
■ sm I 



2 ^'^ ' 



2 1 

sin On, - Z„ 



(AlO) 
(All) 
(A12) 



where g is the SU (2) charge, 9^ is the weak mixing angle, 
and Vud is the CKM matrix element corresponding to u 
and d quark mixing. 

If we define the interactions with background fields as 



(A13) 



define electroweak interactions as 



^ cos6i,„ ^ 



gVud Jli^W+^^ - gVud J^i^W-^" , (A14) 
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and use Eqs. (jA10[) to (IA12I) . we can see that 



tL - 



2 i^ifJ. + ^i/i) ; 



A, 



J, 



J. 



EM 
NC 



T/0 



2 



— sm f „, 



(A15) 

(A16) 

(A17) 
(A18) 



Here, is the baryon current, defined to be coupled 
to v^^^. These relations are consistent with the charge 
algebra Q = T'^ + B/2 (where B is the baryon number). 



V^'^ and A'^ are the isovector vector current and the 
isovector axial- vector current, respectively. We do not 
discuss "seagull" terms of higher order in the couplings 
because they do not enter in our calculations [l^ . 



Appendix B: Form factors for currents 

Here we use matrix elements of the various currents 
to define the form factors produced by the Lagrangian 
Q . By using information presented in Appendix. |X] and 
the Lagrangian in Sec. HIl we can determine the matrix 
elements. 



{N,B\V^\N,A) = 



Uf"f^,u, + Uf{q 7^- ^q^,)u, - _ — u/7 



{B\j\A) 



{N,B\J^\N,A) = 



i)_ CTfii^iq" fp9p q^ _ CTf.^iq" 



2M 



g2_^2"/- 2M 



{B\^\A) , 



Ufjf,u, + -j^ Uf{q 7^- iiqf,)u 



2gv q^9pu - q^qu 



w/7 Ui 



2M 



{N,B■^T,JMA'^N,A) = -tfi {B\'—\A)ufYu 



9p_ q-{q~ K)9tiv - {q- k„)pqi' 

9i 



3^7 q - 
'2fv9v q^ - cTpiyiq" 

9pi^ + -jjj {q-iq- K)gnv -{q- K)pq^) 



{q - - ml 



(Bl) 



(B2) 



2M 

r 



2a(i) - q-jq- fc^) 

g^ {q - KY - ml 



(B3) 



Now we consider {N , B\A^^\N , A) and 
(7V,S;7r, j|V^'|Af,^). In the chiral limit, we find 



{N,B\Al\N,A)=-{B\-\A)un''l 



9 A 9tiu 



qtiqv\ Pa ( 2 \ 
- TH\9 9^1' - q^qi^) 



„ q 9^^ - q^-q^ 

■^Caigai 2 



ml: 



(B4) 



{N,B;7r,j,k,mN,A} {B\-\A}uf^-^j 



'^Ca^gai 



^(1) 

9A9py - {q-iq- K)9py - {q- K)pqu) 



q- {q - K)gf_,y - {q - K)f_,qi, 
{q - k^Y - ml^ 



(B5) 
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Suppose that there is only one manifestly chiral- 
symmetry-breaking term, i.e., the mass term for pi- 
ons; then the pion-pole contribution associated with the 
QA coupling in {N , B\A^^\N , A) will become gA[g^iv — 
q^qu/{q'^-ml)], while the other parts in {N,B\A'^\N, A), 
as well as the whole {N, B; tt, j|y^'|iV, A), will remain un- 
changed. However, we must realize that there are other 
possible chiral-symmetry-breaking terms contributing to 



{N,B\Al\N,A}. For example, (ml/M) Nij^{U - U'')N 
can contribute to {N, B\A'l\N,A) as 



— m': 



{B\^\A) 



To simplify the fitting procedures, we use the following 
form factors [G™'^ can be found in Eq. ([T^ ]: 



{N,B\A^^\N,A) = -GTiq'){B\-\A)uf g,. 



q^-q^ \ ^M„,5 



{N,B-TT,j,K\V;,\N,A) = ^ {B\ — \A)usY^^ 



gAg>.u + SG^X^iiq - K) 



TTu^ , (B6) 

21 g ■ (g - K)g,j,u - {q- K)ii.qv' 



{q - KY 



(B7) 



Finally, we calculate the pion form factor (tt, fc|V?|7r, j): 



(tt, fc, fc,r|V?|7r, i, K- q) ^ ie\{2k^ - q)^ + 2i 



[q ■ Kqp, - q K,,) 



gpiTTV ^ij q_ 1 „ „2; 

57 q^ - 



1 



q- -> TOp in numerator — > ie \ (2fc^ - g)^ + 2i e'-'j, — - (q ■ k^qf^ - g^^Tr^) 



57 g" - 



(B8) 



r 



Appendix C: power counting for diagrams with A 

Including A resonances in calculations, we have a new 
mass scale 5 = m — M w 300 MeV. We must also 
consider the order of the A width F. Formally, it is 
counted as 0{Q^ /M'^); however, numerical calculations 
with Eq. (jD2[) indicate that it should be counted as 
0{IOQ'^ /M^). Because of these two issues, we have to 
rethink the power counting of diagrams involving 6 in 
two energy regimes. One is near the resonance, while the 
other is at lower energies, away from the resonance. In 
the resonance region, the A propagator scales like 



1 



O 



1 

M 
M 



1 



1 



mO{Q^/A'P) 



0{Q)iO{Q) 



1 0{QyM) 
(CI) 



where the 0{1/M) comes from non-pole terms. In the 
lower-energy region. 



S f 



1 



2[5-0{Q)] - 10iO{Q^/AP) 



+ 



1 0{Q) 



1 0(Q) 



0{Q) 26 ' " \M J ' 0{Q) M ' ^^'^^ 

So compared to the normal power counting mentioned 
above, in which the nucleon propagator scales as 1/0{Q), 
for diagrams involving one A in the s channel, we take 
V ^ V — 1 in the resonance regime and v ^ v + 1 away 
from the resonance. 



Appendix D: Renormalized A propagator 

In this work, A's propagator [U is dressed as 
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p2_^2_n(p2)+zTOr(p2) 12 21 3„i2V/' I 22 

+ 0(r/m) X non-pole terms, (Dl) 



2 



r(p') = [p" + M2 + 2A/m) X [(p2 _ m2)2 - (p2 + 3m2)^2] ^(p2 _ M2)2 - Ap^ml . (D2) 



I 

Here: We take m = 1232 MeV as the Breit-Wigner mass [H^ 

and set IT = 0. Note that F is imphcitly associated with 

pd)^"^ = g'"" 7''7'' H J 7[''p''l ^ ^ p^'p" , a factor of e[p2 _ (M + m^)2]. And no singularity exists 

"^P "^P in this propagator at p^ = 0. 

(D3) 

P^P^"" = ^ T^'' - ^ 7 V ^ - ^ P^'' , (D4) 



P^l^^" = ^ (-p^"- + ^) , (D5) 



Appendix E: kinematics 



P^l^^"^ = —p'^p" (D7) Following a standard calculation, we find the total 



22 ~ 2 - - ■ ■ 

P cross section: 



M 



f|2 



M^■Pn^\ 



^ipk-p"" 



li t^ni 



{2^Y2El^^ {2nf2EL i2n)32E^f 

1 d^Pu d^p^ 



2Enf {2TT)^2E[j i2n)32E^ 



|M|2 1 \p,\ \p 



32M„ (2^)5 + 



df}^ dE[f dVtff . 



(El) 



The variables without and "L" superscript are mea- 
sured in the isobaric frame (where A is static). It is quite 
complicated to calculate the boundary of phase space in 
terms of the integration variables in the preceding equa- 
tions. Later, we will work out the boundary of phase 
space in terms of the invariant variables and -M^„ in 
the cm. frame of the whole system, so we would like to 
have the following: 



By using the invariance of the cross section with respect 
to rotations around the incoming lepton direction, we 



have / dTlff = J d cos 9h 2n, and thus 



M 



64M2 {2nfE^+Enf \pj^\E^ 



dn^ dM^,^ dQ 



(E5) 



Ml 



^ Mf J + 2Ek{Etj~\ptf\ cos efj) , 



dQ^dM, 



pL? 



/r2 



2M^{Eh - Eti. 



2 



AM^Eli\ptf\dEf^fdcoseff 



(E2) 

(E3) 
(E4) 



In the isobaric frame, there is no constraint on the di- 
rection of the outgoing pion due to the kinematics. Thus 
the boundary of fJ^ is the whole solid angle in the iso- 
baric frame. Now let's work out the boundary of phase 
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space in the cm. frame. We have 



M 



L\2 



A=PA^ (.Pn^+Pur = (Mn + EfiY ~ {Eli) 



M. 



ML^ip^+pnff = {PA-p'rf) 



C\2 



= MX + Mtf - 2MaE, 



'If 



(E6) 



(E7) 



Here, Ef^j is the final lepton's energy in the cm. frame. 
From now on, all the quantities in the cm. frame will be 
labeled in this way. So, for given Efl, i.e.. Ma, we can 
see that 



M„ + M^ s$ A/^„ i^MA-Mif . 
By using Eq. (|E7p . we find 

Ml + Mff - (M^„)^i„(^ax) 



(-^// )niax(min) 



2Ma 



(E8) 



(E9) 



Then, for given Efl and M^m (or -Ejj), using = 

-Mff + 2EfiE^f - 2Efi\piff cosOff [where Of^ is the 
angle between the outgoing lepton's direction and the 
incoming lepton's direction in the cm. frame, and 
E^^ = {M\ - MI)/2Ma is the initial lepton's energy 
in the cm. frame], we finally arrive at 



[<3^(^^i/)]max 



E^ 



A/5 



, (ElO) 



Ml 



If 
(Elf 



These equations give a description of the phase-space 
boundary in terms of the invariants M-,^n and . 
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